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Abstract 

We discuss the critical behaviour of 2D Ising and q— states Potts models coupled by 
their energy density. We found new tricritical points. The procedure employed is the 
renormalisation approach of the perturbations series around conformal field theories 
representing pure models as already used for disordered spins models. This analysis 
could be useful to understand the physics of coupled critical models like the fully 
frustrated XY model. 
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The critical behaviour of the two dimensional XY-Ising model, consisting of XY and 
Ising models coupled though their energy densities has been the subject of an amount of 
work 0, |]. The model is expected to include in its range of parameters a class of models 
with U (1) x Z2 symmetry as for example the 2D fully frustrated XY model [g], the 2D arrays 
of Josephson junctions [|J, the 2D ANNNXY model || and so on. The phase diagram of 
the XY-Ising model is interesting in its own right because of some unsual critical behaviour 
||. Indeed, there is a plane in the space of parameters with a continuous critical line with 
simultaneous XY and Ising ordering which joins a branch point (where the XY and Ising 
transitions separate) and a tricritical point. Along this critical line, critical exponents asso- 
ciated to Ising-like order parameter were found to differ from pure Ising models. The central 
charge appears to increase continuously from c ~ 1.5 to c ~ 2. Nevertheless, some strong 
finite size effects have been shown and therefore the value c = 1.5 is finally not excluded 
||. Hence, two scenarii are still possible, either the existence of a new universality class or 
simply the superposition of Ising and XY transitions with important finite size effects. 
Whatever the different scenarii, it seems important to understand in a more general context 
coupled crtical models. In a first step in this direction, we wonder in these letter whether 
such puzzling behaviour could also be encountered in systems with Z(2) x Z(N) symmetries 
or more generally with Z(2) x Z(N) n . Namely, could the mixing of two different critical 
models by their density energy create new critical behaviour ? 

When N = 2, the model consists of the n colored Ashkin- Teller model. This model is inte- 
grable when n = 2 |. For n > 2 the transition is either first or second order depending on 
the sign of the energy-energy coupling term [0. In the rest of the letter, we especially focus 
upon the case N = 3 (the case iV = 4 will be a little discussed) because the 3— states Potts 
model belongs to the unitary minimal serie. 

Conformal field theory seems to be the best tool to describe the critical behaviour of such 
2D critical models. Moreover, perturbed conformal field theory has been recently a useful 
tool to investigate the critical properties of the two dimensional Ising and Potts models with 
quenched random bounds || |9], |H| |Tl[]. In this letter, we used the same formalism to inves- 
tigate the critical behaviour of two pure models (Ising and q— states Potts models) where we 
add some energy-energy coupling terms. Such an analysis will also be useful to character- 
ize new simple fixed points, which could serve as tests for pertubation conformal field theory. 



Let us consider the general case of m critical Ising models and n critical Potts models 
which are coupled by their energy density. The Hamiltonian has the following form 

H=J2 H Ln 9 + E H b p otts-Y,9c d j ' d 2 z e c (z)e d (z) , (1) 

o=l b=m+l c,d 
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where c,d belongs to [l..m + n). Hf sing and H b potts represent respectively pure Ising models 
and Potts models at their critical temperatures. e c is the energy operator of an Ising or 
Potts model depending on the value of c. We have a (m + n) x (m + n) matrix of coupling 
constant. We will study the influence of this last term on pure models in perturbation. 



The partition function can hence be written as 

m+n 

m m+n 



- J2 S o,a - Yl S o,b + g c d£ c (z)e d (z)d 2 z . 
Z=U Tr a!S . a J] Tr b , Si b e a = l 6 =-+ 1 ^ d (2) 

a=l b=m+l 

In our notation, a corresponds to Ising models, b to Potts models, and s i c is the spin operator 
of the c th model. A correlation function < 0(0)0(72) >, where O is some local operator, is 
expanded perturbatively like: 

< 0(0)0(72) >=< 0(0)0(72) > + < 5/0(0)00*2) > +^ < tfO(0)O(R) > + • • • 

m m+n 

where <> means the expectation value taken with respect to So = J2 £o,a + J2 S'o.b an d 

a=l b=m+l 

Sj = I Hj(z)d 2 z = J J29cd£c(z)e d (z)d 2 z . (3) 

The calculation of correlation functions can be performed with the Coulomb-gas representa- 
tion fll2"f. The central charge is written as c = h + e, where e will be used as a short distance 



regulator for the integrals involved in correlation functions calculations. In addition, we also 
use an I.R. cut -off r. Then the limit e —>■ corresponds to the Ising model while the Potts 
model is obtained for some finite value of e. The central charge c will be characterized in 
the following by the parameter a 2 + = T^Jry — f + e with 



c = 1 — 24«q ; a± — a ± J aft + 1; (4) 
a + a_ = —1. 

For the pure Ising model a 2 + = | and c = \ while for the 3— states Potts model a\ — |, 
c = g and e = — ^. The 4— states Potts model can be considered as the limiting case of 
this perturbative scheme 0, it has c=l and e = — |. The vertex operators are defined by 
V nm (x) = e m " m *( x ') where $(x) is a free scalar field and the a nm are defined by 

OLnm = -(1 - n)cii_ + -(1 - m)a+ . (5) 
The conformal dimension of the operator V nm (x) is A nm = —(Xjma nm with 

a^r = 2a - a nm = ^(1 + ri)a_ + -(1 + m)a+ . (6) 
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The spin field a can be represented by the vertex operator V p , v -\ whereas corresponds to 
the energy operator e. Note that in the Ising case the a operator could also be represented 
by the 14,1- When we have only one coupling constant go, its renormalisation is determined 
directly by a perturbative computation, g is also given by the O.A. producing 
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go 



£ a {z)e b (z)d 2 z + ^(go j e a {z)e h {z)d 2 z^j H = g J e a {z)e b (z)d 2 z , (7) 



with g = <7q + A 2 g 2 + ' ' ' where A 2 comes from the contraction 

- J e a (z)e b {z)d 2 z J e c (z)e d (z)d 2 z = A 2 J e a (z)e b (z)d 2 z (8) 

This procedure generalizes in a straighforward way to our case. If we suppose that g c d = g\ 
for c,d G [1, to], g c d = g 2 for c,d G [to + 1, m + n] and finally g c d = g\ 2 for c G [l..m], d G 
[to + 1, to + n] or vice versa, we obtain the RG equations associated to these three coupling 
constant (at second order) 

Pl = ^l) = Mm-2) g 2 + 4nng 2 2 +0(g 3 ), 

fa = = - 3eg 2 + An(n - 2) g\ + 4ttto g\ 2 + 0(g 3 ), (9) 

012 = ff 1 ^ = -3^i2 + 4tt(to - 1) g x g X2 + 4yr(n - 1) g 2 g 12 + 0(g 3 ). 

We have introduced a UV regulator e. In order to recover the case of Ising coupled to Potts 
models, we will take at the end the limit e — > — 2= for the 3— states Potts model and e — ► — \ 
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for the 4— states Potts model. Let us study the fixed point structure in the general case. 
When to > l,n > 1 we have only (at this order) two decoupled structures defined by 
gl 2 = 0. It corresponds in the infrared limit to the to— colored Ashkin Teller model || 
plus a n— colored potts model. The critical behaviour of the latter is not properly established: 
for n = 2, the model turns out to be integrable presenting a mass generation indicating a first 



order transition [13] whereas for n > 2 no exact solution exists. Nevertheless a perturbation 
analysis to third order tends to show a first order transition for g 2 > and a new non-trivial 



infrared fixed point for g 2 < ]TT]. At this level, we have not new non trivial fixed points 
mixing both models. 

If we consider now the peculiar case m = n = 1. The coupling constant g\ 2 always flows 
in a strong coupling regime indicating probably a mass generation. Consequently, when 
one superposes a critical Ising model to a critical 3 or 4 states potts model, no new critical 
behaviour is found pertubatively when coupling both model by their energy density. Of 
course, new non-perturbative fixed points cannot be excluded by this method. Therefore, It 
suggests (perturbatively) that no line with simultaneous disordering of the Ising and Potts 
order parameter survives in the phase diagram of this model contrary to those of the FFXY 
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model. Let us now treat the most interesting case m = l,n > 1 as mentioned in the 
introduction. Six fixed points are found: 

9*1 = 9*2 = 9*i2 = (10) 

9*i=9*i2 = ^ 9*2 = 4?r( f_ 2) if n>2 (11) 

* -3erin „ 3e # -3eT 2 ^/n 

Sn{n — 1) 87r(n — 1) 87r(n — 1) 

with Ti,T2 = ±1. There are now four new non trivial fixed point (P"2"|). As usual, we can 
study the stability of each of these fixed point by re-expressing (^]) around the above fixed 
points (<7j = g* + 5gi) and keeping only the smallest order in e. This will give us a linear 
system 5g = A Sg, with 5g a three component vector and A a real 3x3 matrix. The 
stability of these fixed points can be obtained by calculating the eigenvalues of A in each 



of these cases. In this way, it is easy to see that (|10l) is unstable, and that (jTTj) is stable 
in the (#2, #12) plane and marginal in the direction g\. It corresponds to the infrared fixed 
point of n coupled Potts models JTTJ. Concerning the other four fixed points (|12j), two are 
unstable in two directions (when r% = —1) and the order two (ti = +1) are hyperbolic 
fixed points, namely stable in two directions. They constitute new tricritical points. The 

stable plane is defined by the vectors (1, 0, 0); (0, ~^~)^z^P^ \ 1) m t ne (9i,92, 912) space. 
A projection of the flow in the (gi2, 92) plane has been given in Figure 1 for the case n — 3. 
The flow is symetric along the 92 axis. If initially g® > then the flow is always driven in a 
strong coupling regime indicating mass generation as for the simple n-colored Potts model 



TTJ . Nevertheless, if g® < and small, the trajectory will be first attracted by one of the 



tricritical point (depending on the sign of g® 2 ) and then flow toward the stable fixed point of 
the n-colored Potts model. For the case n = 2, the fixed point on the #2 axis does not exist, 
so we always flow in the strong coupling regime except when we are on the stable plane in 
the (<7i, <72> #12) space. It is a reminescence of the exactly integrable 2-colored Potts model 



which is always massive 13 



We want to characterize more accurately these new tricritical points. In this perspective, 
we will compute the renormalisation of the critical exponents associated to the spin-spin 
correlation function. Hence, we need to calculate the multiplicative functions Z ai for the 
Ising spins and Z a2 for Potts spins. A convenient way to define Z a is to add to the general 

n 

Hamiltonian (Q) source terms for the spin operators as h\j ai(z)d 2 z + h\ J <J2(z)d 2 z. 

a=l 

Hence, the multiplicative functions Z a are defined by 



> 2 z 



hxj a 1 {z)d 2 z + h 2 J J2a 2 {z)d 2 {z) = h\Z^l a 1 {z)d 2 z + h° 2 Z a2 J Y,a 2 {z)d 2 

a=l a=l 

As for the coupling constant g^ the functions Z a can be computed in perturbation. The first 
contribution happens to second order. For the details of calculations, we have refered to |§ 
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where rather similar technics have been used. Therefore, the a fields get renormalized as 



02 



a 1 (l + A 2 (^ 2 ) 2 + A 3 (^ 2 )^ 2 u + 



\2 n 



a 2 (l + B 2 (g° 2 ) 2 + CM Vl 



= Z a2 (J2 



with 



A 2 
A 3 
B 2 

C 2 

B 3 

C 3 



-8n(n - l)vr 3 
(n-l)Tr 2 - — 



-6c 



6e 



1 + 2 



r 2 



I)r 2 (|: 



r2 (-D r2 (-ij 



-3c 



7T 



r 2(_|)r 2 (i) 



•127r 3 (n-l)(n-2)- 



-9c 



9e 



1 + 



4 r 2 (-f)r 2 (|) 
3P(-|)r 2 (-i; 



,-6e 



be 



r 2 (-|)r 2 (i 

1 + 2 — ^ 



T 2 (-|)r 2 ( 



(13) 
(14) 



(15) 
(16) 
(17) 

(18) 

(19) 

(20) 



In fact, to obtain such results, we have to conserve two regulators €\ for Ising operators, e 2 = e 
for Potts operators, first to follow the dimensions of expressions and especially because some 
divergent terms (in ^-) appears at one and two loops. These terms disappear when we re- 
place the bare coupling constants by renormalised ones assuring the renormalisability of the 
theory. They have not been included for readability. 

Therefore, using the remormalisation group equations for (^), the equations for Z ai be- 
comes 

■2 



d\og(Z ai 
d log(r ) 

d\og(Z a2 
d\og(r) 



where we have defined 



7T 

-3n# 2 2 (r)_-e + An(n - 1)tt 3 gl 2 (r)g 2 (r) 



7T 



-3(n - l)g 2 n 2 e [1 + 2TZ] - 3gf 2 —e [1 + ATI] 



+4(n - l)(n - 2)ir 3 g 3 2 (r) + 4(n - l)n 3 g'i 2 (r)g 2 (r) 



3„2 



(21) 



(22) 



K 



r 2 (-pr 2 (|) 
r 2 (-|)r 2 (-|; 



The Callan-Symanzik equation gives the form of the correlations functions between spins a. 
For one coupling constant g, we have 



9« 



2 / 



< a(0)a(sL) > 



90" 



e 90 



-2A„ 



< (t(O)o-(L) > r , g(s) 



(23) 
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where we used the notation : ^(rV = laid) > 9( a ) = 9o an d d( s ) * s defined by / (3{g)dg = 

go 

ln(s); r = sa, and a is a lattice cut-off scale. L is an arbitrary scale which can be fixed 
to one lattice spacing a of a true statistical model. The dependence in s of the term < 
a(0)a(L) > r ,g(s) is thus negligible and can be considered as a constant. We want to compute 
the corrections to the critical exponents close at the infrared tricritical points defined in (|12|) . 

So, the integral is dominated by the region g ~ g* and we have / ^^ydg ~ 7s(g*) log(s). 

go (! " 

We thus obtain 

< a(0)a(s) > 90 ~ -r^-^O?*)) (24) 

Therefore, we can compute the corrections to critical exponents associated to spin operators 
dj using the expressions of the tricritical points (|i~2|) 



1<n (9c) = -3n^e(g* 12 ) 2 + 4n(n-l)ir 3 (g* 12 ) 2 g; = + O(e 4 ) (25) 

la 2 {9c) = -3(n-l)(^) 2 7r 2 e[l + 2^]-3(^ 2 ) 2 ye[l+4^] 

+4(n - l)(n - 2)vr 3 (^) 3 + A(n - l)n 3 (g* 12 ) 2 g* 2 
- "27-3 27(2* - ^ 



128(n-l) 2 32(n-i; 
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with e = -r| for the 3-state Potts model and e = for the 4-state Potts model. Hence, we 
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obtain the unavoided result, that there is no renormalisation of the critical exponent A CTl 
associated to the Ising spin at third order in e! It is diificult to give a particular interpreta- 
tion to this result. Yet, the corrections to are non zero and different from those of the 



infrared fixed point of n coupled Potts models [PH . Therefore, these two tricritical points 
are new and non-trivial. We hope that these corrections to the critical exponents of the spin 
correlations functions could be useful to test numerically the results presented in this letter. 
Before concluding, a few remarks need to be done: Suppose we have a theory described by 
one of the tricritical fixed points TJ, i — 1,2 (see Figure 1), then a small perturbation can 
send the flow to the fixed point D characterized by decouple Potts and Ising models. The 
central charge at this point cd is thus the sum of the central charges of each model. Therefore 
according to Zamolodchikov c-theorem [|TJj], the central charge at Ti (ctJ verifies crp i > cd- A 
similar scenario could be imagined for the FFXY model with a new tricritical point inducing 
a strong cross-over regime (it would justify strong finite size effects and cffxy > §)■ 



We have thus shown that the behaviour of coupled minimal models is far from trivial. 
Two new tricritical fixed points have been identified using perturbation around Ising model. 
Maybe, a similar procedure could be investigated for Ising- XY models with another pertur- 
bation scheme. The study of these models with weak disorder is also interesting and wil 
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be presented elsewhere in a longer paper where some details of the present paper could be 
found. 
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FIG.l : The projection of the flow in the (gi 2 , #2) plane for one Ising model coupled to 
three Potts models. Two tricritical points are found. 
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